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Abstract - f < ' ' ‘ ^ ° ' 

'^-This paper is the second in a series of three that analyze a method of 
dimensional reduction. It contains some results for approximation of functions 

4 ) 

on the interval [-1,1] with elements from the nullspace of P , N ^ 1 , where 
P is a second order ordinary differential operator. A special case of this is 
approximation by polynomials. 

The one-dimensional results are used as a tool to prove similar versions in 
several dimensions. These multi-dimensional results are directly related to the 
approximate method of dimensional reduction that was introduced in [13], and they 
lead to statements about the convergence properties of this approach. 

The third paper, which analyzes the adaptive aspects of the method, is 


forthcoming. 










1 . 


Introduction 


In a recent paper, [13], we introduced the concept of dimensionally reduced 
solutions to an elliptic boundary value problem. These are obtained by projecting 
(in the energy) the true solution of the boundary value problem in the n+l-dimen- 
sional domain w x [-h,h] onto spaces of the form 

h N 

V N - { £ W (x)4>.(y/h)| w. arbitrary) , 

J-0 2 2 2 

00 

where {<(>j}j = Q is a given set of functions on [-1,1], (x are coordinates on go 
and y ranges over [-h,h]). For some basic ideas behind this concept, see the 
introduction to [13]. In that paper the focus was on the right selection of the 
(j>j's. It was shown there that for a very wide class of problems the <J>j's should 
be selected such that 


spanUj}^ 1 * N(P k ) , 

where P is a second order differential operator intrinsic to the elliptic boundary 
value problem. 

The estimates of the error given in [13] were asymptotic in h ■* 0 . The present 
paper, which was already announced there, treats convergence as N -*■ ® for a fixed 
value of h . For convenience the fixed value of h is set equal to 1. 

If the bilinear form associated with the elliptic boundary value problem 
satisfies some kind of "inf-sup" condition, then it is well known that the rate 
of convergence is the same as the rate of approximation (cf.[l]). 

The results proven here are hence formulated as approximation-theoretic 
estimates, and as such have Interest regardless of the concept of dimensionally 


reduced solutions. 











2 1 

The results are all concerning approximation in the L - and H -norms, 

i.e. ideally suited for second order problems. This is not crucial and similar 

results can also be obtained, e.g., for the norms (/ 1 ||4^|| 2 dy + J 1 !|a 1 ^ 2 u|^ dy) 1 ^ 2 

-1 dy -1 

introduced in [13]. (A here denotes a strictly positive definite (unbounded) 
operator in a Hilbert space H , and u is a function with values in H .) 

For reasons of convenience the approximation results are formulated without 
any boundary conditions. Various types of fixed boundary conditions can immediate¬ 
ly be included based on the present proofs. 

Estimates of the error introduced by dimensional reduction, as N goes to », 
do exist in the literature (cf. [5,7]). The problems considered in those two 
papers come from structural mechanics. The elliptic operators have constant 
coefficients, i.e. the <j>j's are polynomials. The results are not nearly as 
strong as the ones established here. In [7] the estimates are based on the degree 
of regularity in C -spaces; this is not very well suited to the regularity pro¬ 
perties of solutions to elliptic boundary value problems, and therefore gives 
crude estimates. The estimates in [5] are based on bounding the remainder in the 
N-th order Taylor expansion. The estimates are very crude and do not give any 
indication of the rate of approximation. 

We now give a short review of the contents of this paper. In section 2 it 
00 

is shown that the set U N(P k ) (N denotes the nullspace) is dense in H^" for 

k=1 d d 

any second order operator P » b a , where both a and b are bounded from 
above and away from 0 . This is the obvious generalization of the fact that the 

polynomials are dense in , and it also justifies the claim that the dimensionally 
reduced solutions introduced in [13] will get arbitrarily close to the true solu¬ 
tion. In section 3 the rate of approximation using functions in N(P^), N >_ 1 , 
is linked to the regularity of u in spaces of the type P(p m ) . This general 
result though is not always optimal, as shown, e.g., by Theorem 4.1. Section 4 
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is devoted to giving a necessary and sufficient condition for a certain rate of 
approximation by polynomials (i.e., the case where the operator P is a constant- 
coefficient operator). In section 5 this is carried over to results in several 
dimensions — directly relating to the concept of dimensional reduction. The exam¬ 
ple treated in section 6 is of the same type as the numerical examples in [13]. 
Finally the appendix contains the proofs of several results about the eigenvalues 
and eigenfunctions for two-point boundary value problems, as used in sections 2 
and 3. 

Note: Unless otherwise stated, all constants denoted by capital letters 
are generic. 
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2. A Density Result 

Let a and b be two functions in L°°([-l,l]) such that 3 constants 

a , b with 
o o 


0 < a Q <_ a(y) 

0 < b 0 1 h(y) 

By P we denote the differential operator 


, d d 
dy a dy 

P is considered as a mapping L^([-l,l]) 3 P(P) -*• L^([-l,l]) . N(P k ) denotes 

the null-space of the operator P for any integer k 1 . It is easily seen 
k 1 

that W(P ) C H ([-1,1]) . The first theorem in this section proves the density 
of a certain class of functions associated with the operator P . 


Theorem 2.1 


0 N(P k ) is dense in ^([-1,1]) 
k»l 


Proof 

- Y [ 1 

By a change of variables, y' = I 

-1 


operator P transforms into 


and multiplication by -1 the 


d £ d_ 
dy' b dy' 


We can therefore for the proof of this theorem assume that P 

- a(y) 4 - where a satisfies: 3 a constant a with 0 < 
dy dy o 

Define the operator Q by 

V(Q) - V(P) fl H 1 ([-l,l3) and Q-P on P(Q) . Let f Q 


is given by 

a 0 la(y) • 


denote the function 
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£ 0 W 


ds € N(p) 

-1 


and define the sequence 


, , 00 

f i i*0 


by 


f ± = Q _i f 0 € W(P i+1 ) . 


0 < \ o < ^ < . . X^ X ^ . . denotes the eigenvalues of Q (repeated 

OO 

according to multiplicity). Let ^ u m ^ m =o an or tbonormal basis of eigenfunctions, 
u^ corresponding to X^. f^ can then be expanded as 


f 


0 


00 


l 


a u 
m m 


» 


and with this notation 


f 


i 


m =*0 


a X -i 
m m 


u 

m 


We now proceed to prove that any eigenfunction u^ can be approximated from 

within (j N(P ) . The proof is by induction in m , and we start with m * 0 . 
k-1 

For any i > 1 we have that 


I l u 0 * “o 1 x o f iN 2 i 1 

H 

C||Q 1/2 < Uo - *J*1>II*2 


/ \2 

/ \2i-l / 

(“jN 

(v*}) - c (v4) 


k2i-l 


l * 
1-1 J 
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where we have used Lemma A.3 to guarantee that a. ^ 0 . Since £ < » 

J-l j 

and, by Lemma A.l, X^/X^ < ^ tk is shows that 


- 1 . 1 - ^ 
a o x o f i " u 0 

u £ u W <P k ) 

k-1 


as i -*■ « , in ^([-1,1]) , 


- denotes the closure in H ) . 

Now assume it has been proven for some m > 1 that 


{u } m_1 C U N(p k ) 
1 j j-0 ~ k-1 ( } 


From Lemma A.3 we know that a ^ 0 , and hence for any i > 1 

m — 


00 . 

u-a 1 X i f = x ~y / a /a \(\ /X V u. 
m m m i m,i ^ j m jj j 


where x 


- - y (a /a, \( X /X A u € U N(P^ 

m.i j-0 V J m J[ m J/ J k=1 


due to the induction hypothesis. As before the H - norm of the sum 




can be estimated by 


(^W“ (L •;)' 


Because of the facts that 


(,L -It 


j-m+1 


< « and, by Lemma A.l, X /X ._ < 1 

m m+i 


this shows that 




•SJW-* - 


“M* 

x.+aXf. -*■ u as 
m, i m m i m 


i -V ao f 


in H a ([- 1,1]) , i.e. 


{u }“ 0 £ 0 W(P k ) 

3 3 k-1 


This finishes the induction proof, and we conclude that 


{u }" C u W(P k ) 
3 3 k=l 


From the definition of Q it immediately follows that 


0(Q 1/2 ) = H 1 (f—1,13) 


- - GO 1/0 

and since iu^J j = q is complete in "(Q ) this proves that 


H^t-l.l]) £ u w ( pk ) 
k-1 

Now if u€H^([-l,l]) we shall, by choosing c - u(-l) and d = 
1 -1 

(u(l)-u(-l))// a 1 (s)ds , obtain that 
-1 


u - c - df Q € H ([-1,1]) 

Since 1, f^ € N(P) we see by a combination of this and the previously 
proven inclusion that 


^([-1,1]) - 0 N(P k ) 
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Based on Theorem 2.1 we can easily prove a result concerning the dimensionally 
reduced solutions as introduced in [13]. This result guarantees the fulfilment 
of the goal stating that the dimensionally reduced solutions shall be able to get 
arbitrarily close to the true solution. 

Let u denote a domain in IR n with a Lipschitz boundary. 

Theorem 2.2 

The set 


W -I i 

{ £ v. (x)q. (y) | J € IN,v, 6H Coo) and q, € U W(P ) for 0 £ j £ J } 

j=0 2 2 2 2 k-1 

is dense in H^(w x [-1,1]) . 

Proof 

Follows immediately from Theorem 2.1 and the fact that 

J , , 

t ^ v. (x)w, (y) | J € lN,v, €H (w) and w.€H ([-1,1]) for 0 <_ j £ J } 
j»0 ^ 2 2 2 

is dense in x [-1,1]) . 

a 
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3. Estimates of the Rate of Approximation 

In the previous section we proved the density of a certain class of functions 

In this section we shall prove 
some results concerning the rate of approximation. The first theorem is the 
following. 


associated with the operator P = b ~ a . 


Theorem 3.1 

2 

Assume that a,b€C ([-1,1]) , and let m be an integer 0 . For any 

e > 0 there exists a constant C such that 

£ 


inf 

v€N(P N ) 



V(? m ) 


Vn > 1 


Note. ||*|] denotes the norm ||P™(*)|| 0 + | | *| | „ 

- P(P m ) L 2 L 2 

One can of course combine the statement in Theorem 3.1 with interpolation by the 

K-method (cf. [4]). This way it follows, that if u 6(L 2 ;P(P m )) for some 

s 

0 < s < 1 , then for any e > 0 


inf 


u-v 


N. 


0 < C N 
. 2 — e 


-ms+e 


v € W(P ) 


(L 2 ,0(P m )) 


s,°° 


2 

The smoothness requirement that a,b€C ([-1,1]) is not necessary; as it 

2 

immediately will follow from this proof we only need that a/b is a C -function. 
This last remark applies to all of the results in this section. 


In order to prove Theorem 3.1 we need an auxiliary result concerning uniform 
approximation by polynomials. This result can be found, e.g., in chapter 6 of [6]. 
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Lemma 3.1 

Let 4> be a function in C^([c,d]). Define $ by $(t) = 4>(—^ cost + , 


t € [0,ir] • 


Let r be a non-negative integer. There exists a constant such that for 


any <f> with $€C r ([0,ir]) the following estimate holds 


inf|*-P N l 0 < C r (N+l)“ r |$| r VN > 0 

Px, 


The infinum here is taken over all polynomials P N of degree £ N . |•| 

i 0 r 

and |*j r denote the norms in C ([c,d]) and C ([o,tt]) respectively. 

We now continue with 
Proof of Theorem 3.1 

Like in the proof of Theorem 2.1 we may also here assume that P is given by 

P =- ^ a(y) . Let f i i > 0 be defined as in that same proof. For any set 

N 

of coefficients {c^}^_q we ^ ave » using Lemma A. 3, that 


^ a 1 
u ~ 1 c.f = l a (8,/a - l c.X'Su. , 

i*0 1*0 3 3 3 i-0 3 3 


00 

where £ B.u, *- s the expansion corresponding to u . If by p 

j=»0 3 3 N 

J r i 

polynomial p (x) = £ c.x * then the above can be rewritten as 

N i=0 1 


If by p„ we denote the 
N 


" ' J 0 c i f i ' £ W a J ' P » ( ^ 1))U J 


and this leads to the equality 


i u - I c i f iii 2 2“ i a i (e i /a i 

i-0 1 L j-0 3 3 3 


P N (A j 1))2 
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As In the proof of Theorem 2.1 Q denotes the restriction of P to 

P(P) D H^([-l, 1]) . Let us now for a while assume that u€P(Q^) . 

— 2 00 

Choose A so that {A^ }^gC[0,A] . Define a sequence of functions 
$ m €C°°([0,.A]) , 1 <_ M, with the following properties 

“ S j /0t j ° - j - M_1 
*=0 on [O.A^] 

Let $ denote the mapping 

$(t) ■ y (1-cost): [0 ,it ] -*> [0, A] , 

it then follows from Lemma A.l and Lemma A.2 that 

|* -1 (A"^ 1 ) - <t' 1 (A" 1 )| > C/j 2 for any j > 1 . 

This estimate tells us that it is possible to construct the <(> M 's such that 

Vm, j >_ l 


IV^CO)! < C sup |(j+l) 2r g,/a 
0<j<M-l J J 


Now since u€P(Q k ) we know that 16j | £ C^Cj+1) 2 ^| |Q k u| | 2 , 

L 

this with Lemma A.4 we get 


Bj/Ojl ±C k (j+l) 


-2k+l, ,„k. 


Q u 


and combining 










i.e., we have for any r >_ k and M ^ 1 


i V $(t)) l r 1 C r M 2(r " k)+1 |!Q k u|| 2 

L 

Because of Lemma 3.1 we can now, for any r ^ k , M ^ 1 , find 
of degree <_ N such that 


IV P n'o 1 C r (N+l)' r M 2(r “ k)+1 j|Q k u|! 2 

L 

We now go back to estimate the right hand side of (1) 

GO M-l oo 

I a 2 (e /a -pjju" 1 )) 2 < l a 2 ((<^ -p^CXT 1 )) 2 + 2 £ & 2 + 

j=0 J J J w 3 j=0 3 w 3 j=M 3 

OO 

+ 2 l a‘((* m-pShxT 1 )) 2 

j=M 3 3 

The first and the third sum can be estimated by 

C r (N+l)- 2r M 4(r_k)+2 ^ a 2 ^||Q k u|| 2 2 < 

< C (N+l)" 2r M 4(r " k)+2 ||Q k u|| 2 , 

L 

the second by 


C 


k 


V (j+D" 4k ! lQ k u| i 2 2 1 C.M‘ 4k+1 
j-M L 



polynomial 
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In summary we have therefore proven 


N 


11“ “ I c f || < C (<N+l)“^ r M' 

i-0 L 


- 2r„4(r-k)+2 + M ~4k+1^|jpk^ |2 

L 2 


for any r ^ k , N ^ 0 and M ^ 1 . Taking M * + 1 ([• ] denotes the 

integer part) this estimate gives 


( 2 ) 


1“ - l c i f il| 2 2 -- C k ( N +1) " 2k+1 |l Qk u|| 2 2 . 

i=0 1 T K * z 


L 


k 2 

all provided that u€V(Q ) . Let denote the L -projection onto linear 

combinations of the functions f^..., f N . 

(2) expresses that 


M n u| | 2 < C k (N+l)' k+1/2 l |Q k u| I 2 , 

L Li 


at the same time it is clear that 


1 l u " n N u l I 2 - I l u N 2 

L L 

Applying interpolation by the K-method we get for any 0 £ m k 


-(k-l/2)f 

il u -n N u|I 2 < C k ( N+l) (|Q u| I 2 

L L 

Now let m be fixed and k -► °® , from the previous inequality we then get 


Ve > 0 3 C 


such that 
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Based on Theorem 3.1 we can prove the following result concerning approximation in 
. 1 

the H -norm. 


Theorem 3.2 

2 

Assume that a,b€c ([-1,1]) , and let m be an integer 1 . For any 

e > 0 there exists a constant C such that 

e 


inf 

v€N(P N ) 


u-v 


H 


, < c 
1 — £ 


N 


,-mfl+e 


ll«l 


V(P m ) 


Vn > l 


Proof 

*\j N 

From Theorem 3.1 it follows that there exist v^€N(? ) such that 



-nH-l+e 

N 



Now choose v fN(P^) with Pv„ * v„ , and such that v„ = u for y * + 1 . 

N N N N — 

It then follows that 


lu-v v 


j < C||Pu-Pv N | 


cliPu-v. 


N' 


0 < C 
2 - e 


N 


-m+l+£ 


u 


P(P m ) 


□ 

We can also easily prove a result relating to the dimensionally reduced solu¬ 
tions as introduced in [13]. Let wciR n be a domain with a Lipschitz boundary, 
x * (x^,...,x ) denotes coordinates in w and y ranges over [-1,1]. P denotes 
b(y) -jp a(y) jp considered as an operator L^ x f-l.ll )3f7(fr) ->■ L^(u x [-1,1]) . 
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Theorem 3.3 

2 

Assume that a,b€C ([-1,1]) and let m be an integer >_ 1 . Let u be 

an element of L (w x [-1,1]) with u,..., y u , u€P(P m ) . Then for any 

X 1 n 

e > 0 there exist C £ (independent of u ) such that 


inf 

v€V 


u-v 


N 


H (m x [-1,1]) 


< r°* 1+ i I n£ 


l L \ 

\i=l 





Here V, 


N 

' denotes the set {£ w ,(x)<t> (y) |w, CH 1 ^)}, where {, J? 6 " 1 is a basis 

j=0 J J J J J =u 


for N(P k ) C. H 1 ([—1,1]) . 


Proof 


Let v N denote the 


N 

orthogonal projection of Pu onto ( [ w,(x)<t >, (y) |w €L 2 (w) } 

j=0 J J J 
3v„ 


/ M A 

in the L (ui x [-1,1]) inner product. Then it is clear that ^ is the L 

3 *i 

projection of ll u onto the same subspace. From Theorem 3.1 we immediately get 


l I I* ■dr (u_v N> M 2 

i-1 dX i W L Z (mx [-1,1]) 


+ Il^(u-vg)I| 2 1 

l/(«x [-1,1]) 


< C^N+l)-^ 1 ^ 



£>(?“) 



for any function with Pv* = v N (if N is odd such a v* will be 

contained in Vj^ * but this is not necessarily so for N even). Now choosing 
v j$ so that also v* ■ u for y * + 1 (this is obviously possible), it follows 


that 
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4. The Constant Coefficient Case - 1 Dimensional Results 

The following two sections are devoted to the case where the operator P 

has constant coefficients. In the previous section we estimated the rate of 

approximation for general P's , but the estimates established there do not have 

exact inverse counterparts nor are they always optimal. As will be shown in this 

section and the next, the question of approximation rate can be much further 

clarified when P is a constant coefficient operator. We start with an analysis 

of the 1 dimensional problem. The space N(? ) , k 1 , consists simply of all 

polynomials of degree £ 2k-l . Theorem 3.1 combined with interpolation says 

that if u€H t ([-l,l]) , then there exist polynomials p^ , of degree N , such 

that I Iu—p„I I _ < C (N+l) t ^ +£ . Under the present simplified circumstances 
N 1 ' z — e 

la 

we can prove a better result. In the formulation of this result we use the Besov 
spaces B~ , t > 0 , (cf. [4]) , instead of the ordinary Sobolev spaces H t . 
For an interpretation in terms of the spaces H 1 " use the inclusions H C C bJ: c_ 
H t-£ valid for any t > 0 , e > 0 . 

Theorem 4.1 

Let t be a given positive number. There exists a constant 
for any uCB^j 0# ([—1,1]) one can find a sequence of polynomials 
degree of p N <_ N , with 

'Kll 2 i c t (iffl) ' t ii u il t 
L B 2,® 

Note. A similar theorem is also valid for the H^-norm. The estimate here 


C t such that 
^ P N } N-0 ’ the 


u-p 


N 1 


, 1 ± 


c t (iw-i) 


-t+1 




' 2 ,‘ 


becomes (for t > 1) 










The rate of approximation established in Theorem 4.1 is optimal in the 
following sense. 

Theorem 4,2 
2 

If u€L ([-1,1]) and there exist a constant C and a sequence of poly- 

ao 

nomials ^ e 8 ree P N <_ N , such that 

||u-P N || 2 £ C(N+1) C for some t > 0 , 

L 

then h€«5,J loc nB|« . 

Note . Theorem 4.2 is not an exact inverse of Theorem 4.1 since it only 
t/2 

guarantees that u€B_ ([-1,1]) . But based on Theorem 4.2 we conclude that 

z,°® 

for a general type function in b!: ([-1,1]) we cannot expect more than an approxi 

mation rate of (N+l) t . 

Theorem 4.2 is optimal in that one can find u such that | |u— p n | | £ C £ (IH\L) 
and u^B^ e , u ^ ( B 2 + «)]_ oc for an y e > 0 (cf. [12]) . 

The proof of Theorem 4.1 is very simple, based on transforming u into a 
periodic function and estimating the remainder of the k'th order Fourier expansion 
Details can be found, e.g., in [2]. 

The proof of Theorem 4.2 is not quite as simple. The cornerstone is the 
so-called Bernsteins inequality 


^dy^N I ^2 



valid for any polynomial of degree <_ N . For more details see [2] or [10]. 

As already noted Theorem 4.2, although optimal, is not an exact inverse of 
Theorem 4.1. This can be taken as evidence that the standard Sobolev or Besov 


spaces are not very good for expressing the kind of regularity needed for a 






certain rate of approximation by polynomials. They do not take into account the 
well known fact, already noted by Timan (cf. [9]), that the polynomials have a 


certain ability to absorb singularities at the end points of an interval. 


d 2d 

Let L denote the operator - -jj— (1-y ) with domain of definition 


P(L) - {u£L 2 ([-1,1]) i Lu€L 2 ([-1,1])} 


Now introduce the Besov spaces H , t > 0 , by 


h 1 * a?a p );m q )) 

where p,q are two integers with 0 <_ p < t < q and 0 < s < 1 is selected so 
that p(l-s) + qs » t . Because of Theorem 14.1 in [8], which says that 

„ - V(L^^ e ) + £ 0 ) > and the reiteration theorem on p.50 of [4] 
it follows that modulo equivalent norms H t is independent of the choice of p 
and q . 

We are now in a position to characterize completely the regularity needed 
for a certain order of approximation by polynomials. 

Theorem 4.3 

Let t be a positive number. For any u€H*" we can find a sequence of 
00 

polynomials ^Pjj^n* 0 ’ t * le degree P N ± N , such that 

I I U “P N I I 2 - l u l I t 

L H 

2 

On the other hand if u€L ([-1,1]) and there exists a constant and a 

oo 

sequence of polynomials {Pfj^N*o * ^ e 8 ree Pfj — ^ * suc ^ that 
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Il u -P N l! 2 1 c u < N+1 >“ 2t 
L 

then u€H t and 


I Ml t £C(C u +||u|| 2 ) 

H L 

for some constant C independent of u . 

Note . C y here is not generic, it is the same constant in the two inequalities. 
Proof 

We start by proving the direct part. It is well known that the eigenfunctions 
of L are the Legendre polynomials {JL • Also 


L(^) - kOefl)^ 


Let now u be an element of V(L P ) and let £ a £ be the Legendre series 

m*0 m m 


for u . Since u€V(L P ) we know that 7 a^m^ < ||u|| 2 

_r. m n 


N 

l « £ » then 

m m 

m»0 


m”0 


P(L p ) 


Define p 


N 


|u-P H l 


m*N+l 


a 2 < 
m — 


(N+l) 


-4p 


m-N+1 


2 4p 
am 
m 


< (N+l) 


-4pi 


I 2 

V(L p ) 


i.e., I Iu-p„| | - < (I«-l)" 2p | |uj | 

N L 2 P(L p ) 

Interpolation applied to this gives the desired result. 

We now turn to the inverse. Assume that there exists polynomials p N of 
degree £ N such that 


r 
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> U_P nI 1.2 - C u (N+i) 

L 


-2t 


Define K - P, , K ■ P -P , m > 1 . 
o I’m 2 m 2 m- ^- — 


Then 


I I* I I 2 - c u + I l u l I 2 and 

L L 


I Kll , 1 | |u-p _| 1 , + | |u-P „ < C/C2 m > 1 


,-2tm 


“ M L 2- "“W 


2 m-l‘'2 - °t u 


2 

Since || Lp || , < Cn || p | | for any polynomial p of degree £ n , it follows 
n L n L n 

from above that for any non-negative integer q 


I |i 9 * 0 l I 2 - C q (C u + INI 2 > and 

Xj l 


[ | || , < C -C.-C .2 2 ^ _t)in m > 1 

11 m' ' 2 — q t u — 


Now define v, = V fi . We then get 

k n m Q 

m=0 


|v ‘WJo <IK M 2+U *-V 


1 C (C u+ Hull + C U l 2 2 ^*) 
T m*l 


< C 2 2(q_t)k (C. + Hull „) 


u " "t 2 


provided q > t . At the same time 







- 


l u -v k U L 2 = 1 lu ~ P 2k l > l 2 -C u 2 


By defining = 2 q we therefore 


(||u-v k || + St ||v k || -)<C (c„4. Hull ) , 

L o 1 g \ ’ L 


p(L q y 


and since s, -*■ 0 for k ® this proves that 


u €(L 2 ;P(i- q )) t - 


Ml 1 1 C cc + ||u[| 2 ) 

H t q,t u 


This theorem also allows a version formulated by using the spaces V(L*~) 
instead of the corresponding Besov spaces. It is derived from the inclusions 
0(L C ) C_H C C 0(L C e ) valid for any t > 0 , e > 0 . 

A theorem similar to Theorem 4.3 but concerning approximation in the 
norm can be derived immediately based on Theorem 4.3. 

For practical purposes, in determining the rate of approximation, the 
following characterization of P(L q ) (cf. [3j) will often be convenient: 


P(L q ) - {u€L 2 ([-l,l])!u€H q ([-l,l]) , (l-y 2 ) q u€H 2q ([-l,l]>} 
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Let us now give one simple example that shows how a result similar to 
Theorem 4.3 can be established also in a case with non-constant coefficient. 


Example 4.1 

Let P denote the operator a 1 a ^ » where the function a is 
given by 


a(y) = 



for y ^ 0 
for y < 0 


with a + and a^ being two positive constants. (This is the operator arising 
in the numerical examples of [13].) 

It is not difficult to see that the following set of functions is a basis 
for N(P N ): 


K 2k 


♦o " 1 ’ *i = 1(71 y 

! f 7 

^2k-l (t)dt ’ ^2k+l a (y) ^ 2 k (t)dt 

-1 -1 


1 < k < N-l 


Here denotes the Legendre polynomial or order k . Performing the 

Gram-Schmidt orthogonalization on the set <f> Q , <fr lt ... ,<j> 2N _ 2 **2N-1 ^ in that 
sequence), in the inner-product <u,v> = u(y)v(y)a(y)dy , we end up with 

a new set of functions i|)^,...,i(> 2 2 * ^2N-1 ' ^k ds a P* ecew i se polynomial 

of degree k . Let denote the operator -a x (y) a(y)(l-y ) , it is 

then clear that 




a y 2k 


i 


k 

l k 
=0 2 


"2J 


0 < k < N-l 
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^a^2k-l * ^ A j,k ^2j-l 1 - k - N 


Now we have because of the orthogonality of the and the fact that 


L a is self-adjoint in L ([-1,1],a(y)dy) 


<L a\’' P ]\ = < \’ L a^j > a = 0 for J < k 


i.e., 


^a^k = W for any k 


It immediately follows that = k(k+l) . 

oo 2 co 

It also follows, since is dense in L ([-1,1]) , that ^k’\^k=0 

is a complete set of eigenvalues and eigenfunctions for L . 

As In the proof of Theorem 4.3 we now get that 


inf 

v€N(P N ) 



< CN 


-2t 


u €(V(L p J; VU q )) a 

a a w i 


for any 0 <_ p < t < q , and 0 < s < 1 chosen such that t = p(l-s) + qs . 

In summary, we have found a singular operator L that characterizes the rate 

a 

N 

of approximation with functions in W(P ) the same way that the Legendre 
operator does with polynomials. 
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5. The Constant Coefficient Case - Dimensions Higher Than 1 

In this section we prove a result relating to the dimensionally reduced 

solutions introduced in [13]. We give a characterization of the regularity 

needed for a certain rate of approximation. The approximating functions are 
N 1 

of the form £ w -«(x)P.i(y) where w. €H (tu) and p is a polynomial of 

j»>0 J J 3 J 

degree j , i.e., the operator P has constant coefficients. w as before 
denotes a domain in IR n with a Lipschitz boundary and y ranges over [-1,1]. 

In the proof of the main result in this section the following lemma will 
be very useful. 

Lemma 5.1 

Let C_ Hq be two Banach spaces with norms |[•[and ||*|Iq respec¬ 
tively. 

Let w =0 be an increasing sequence of subspaces of and let 6 be 

a positive number. We assume that the following implication holds 

u€H and inf | |u-q| L < C (N+l)” 6 , Vn > 0 

q€V N 

d- 

u€H l and M u l l x 1 C(C u + | |u| | Q ) 

for some C independent of u. (C^ here is not generic, it is the same con¬ 
stant in the two inequalities.) 

As a result of this it follows that for any 0 < 0 < 1 

u€H q and inf ||u-q|| 0 <C (N+l)“ 96 , VN > 0 
q€V N 

u€(H o ;H i ) e,» and IMI 9 ,»i c C c u + IMI 0 > 
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for some C independent of u . (As before is not generic.) 

Proof 

Let 0 < 0 < 1 and assume that there exists a sequence of elements 
q N^ V N’ N ^ 0 , such that 


Define 


I l u -q N l l 0 i c u (n+d 


-06 


*0 = q l 


K = q -q , m>l , 

m 2 m ^m”! — 


then 


I" - l *.llo ■ i c „ 2 

m»0 2 


-k0B 


At the same time 


1*01 <0 - C u + I '0 and 


Mi 0 i h^ 2 ji 0 + n -,^11 < c 6iS c u 2 - 86 , „ > i 


That is, h. €V m > 1 , and 

m 2 m “ 


2 


” S( 6 ' 1 \lloiCe, 8 «u + H“Ho > 2 " 6 


From the first implication in the statement of this theorem it follows that 


|2* s(0 - 1, \ 1 II 1 «c 9>6 (C u + ||u|| 0 ♦ ||2“ 6 < 9 - 1 \|| 0 ) 



Ml i c e,B <C u + ||u|| 0 )-2 m6<1 - e) 


We therefore get 


kB(l-6) 




If we define =■ 2 the following inequality now holds 


S k (||u " J 0 ! 0 + Ml> - C 0,B (C u + N«ll 0 > 


Since s, -*■ 0 for k -*• ® this proves that 


u€(H () ;H 1 ) 0>oo with 


l«ll e ,.iC et6 (C u + ||u|l 0 ) 


Let us introduce the spaces 


K* « {u€0(L R )|-^u€P(L R ) i-l,...,n} , 

3 2 3 *y 

L here denotes - (1-y ) considered as an operator L (w x [-1,1])3 

2 

V(L) ■+■ L (oj x [-1,1]) , and R is a non-negative integer. 

N 

V N denotes the space { £ w (x)p (y) |w. €H J '(u>)} where p is a poly- 

j«0 J J J J 

nomial of degree j , j 0 . We are then able to give the following character¬ 


ization of approximation by the spaces V„ in the H-norm. 

N 


i 





Theorem 5.1 


1 R 

Let a be a given positive number. If uf(H (u x [-1,1]); r) a/R,» for 
some integer R > a and R > 2a/e , where e is a positive number, then there 
exists a constant C such that 

inf [|u-qjj < CCm) -2a+£ , VN > 0 . 

q€V N ir(u)x [-1,1]) 

On the other hand if for some e > 0 there exist a constant C such that 


then 


inf | | u-q j | 1 C(N+1) 2 “~ e , \/N >_ 0 , 

qCV N H ( ci) x [-1,1]) 


u6 n 

R€IN ,R>a 


(H 1 (o) x [-1,1]) ;K R ) 


a/R,‘ 


Before we proceed with the proof of Theorem 5.1, let us state a corollary 
that immediately follows from this theorem. 

Modulo e this is the equivalent of Theorem 4.3 in more than one dimen¬ 
sion. 


Corollary 5.1 

Let a be a given positive integer. If u€ fl (H 1 (cd x [-1,1]) ;K R ) . 

R€lN,R>a a/K ’ 

then for any e > 0 there exist a constant C £ such that 


inf | | u-q | | , _< C (N+l)~ 2 “ +e , N> 0 . 

q€V N H^Co) x [-1,1]) 


On the other hand if for some e > 0 there exists a constant C such that 





-2a-e 


inf l|u-q|| H l (ajx { _ 1M < C(N+1) , V«iO 


u€ n (H 1 ^ X [- 1 , 1 ]);K R ) 
R€IN,R>a 


Proof of Theorem 5.1 

Assume that u€K and R is an integer > 1 . Let v denote the 

^ N N 

orthogonal projection of Lu onto { £ w (x)p . (y) |w. € L^(w) } in the L z (cu x 

^ =0 3 2 

[-1,1]) inner product. Then it is clear that v N is the L projection 

of L - 5 — u onto the same subspace. From Theorem 4.3 we immediately get 


n ^ 


J, | I L ^T < “- V N )II i 2. , . + l|Uu - V N ) l' I 2. . . ... 

1=1 i L (a) x [-1,1]) L (to x [-1,1]) 


< C r N 


^ + Hull \ 

i -1 9x i V(L R ) P(l R ) 


for any function v*€V^ with Lv* » v^. Now choosing v* so that also / v*(x_,y)dy 

,1 _1 

B J u(x,y)dy for any jc€w (this is obviously possible) it follows that 


Jl^3x (u-v Kpl I 2. , n + ^ U-V N^ 2, r n + I l3y (u_V N^‘ 2. r i its ~ 

i*l i L (a) x [-1,1]) L (w x [-1,1]) 7 L (oj x [-1,1]) 




+ I |Uu-v*) I I 2 ) 1 

1 ]) N L^w x [- 1 , 1 ])/ 


U| I + I|u|| -w 

V(L ) V(LT' 


,) ■ c * ” 2 


i 






31 


Using interpolation on this result we get that 


2 a 

inf | | u-q | | < C N _2ct+ R | 

?€V n H (oj x [-1,1]) 



where lll u lll a / R a> ^ enotes the norm on (H^(a> x [—1,11 ;K R > a /R • Since 

2a/e <R , i.e. 2a/R < e , the direct part of this theorem immediately follows. 

Let us now give a proof of the second part of the theorem. If 

-2R- Re 

inf ||u-q ]| . £ CN ""®" for some R > a , then as in 

q£V N H X (u) x [-1,1]) 

the proof of Theorem 4.3 it easily follows that 


u €V(1 R ) and u €V(L R ) , 

i.e. u€K R . If we apply Lemma 5.1 with * K R , - H^(u x [-1,1]) and 

9 =* a/R , we then get that 


inf ||u-q || < C(N+l)' 2a ' E 

q€V N H X (<i> x [-1,1]) 


implies 

u€(H 1 (o) x [-l,l]);K R ) a/R>00 , 


for any integer R > a , i.e. 


u € D (H 1 (a> x [-1,1]) ;K R ) „ 
R€DJ,R>a 0/ ’ 


□ 


i k 
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For the conclusion of this section let us give a simple example that shows 
the practical usefulness of Theorem 5.1 (or corollary 5.1). 

Example 5.1 

Let w be the interval [0,1] . Let y be a positive number and let 
(A.,9) denote polar coordinates around the point (1,1). We then consider func¬ 
tions of the type u ■ A Y $(9) , where <t> is an element of C°°([0,ir/2]) . 

It is not difficult to prove that 

u€ 0 (H 1 (to x [-1,1]);**) . 

RGN ,R>t t/R ’“ 

for 0 < t <_ y , and that for a general choice of 4 > this is not so for any 

t > y (if Y ^IN, then this is not so for any <f> and x > y except = 0). By 

an application of corollary 5.1 we therefore get that 

inf | |u-q| | < C (N+l) -2Y+e , VN > 0 , 

q£V N h(ui x [- 1 , 1 ]) E 

for any e > 0 , and at the same time that for a general choice of <J> (or for 

any d> ^ 0 in the case y^IN) there exist no t >0 and C £ such that 

inf I |u-q j | < C (N+l)" 2Y_e , VN > 0 . 

q€V N H (w x [- 1 , 1 ]) 

y 

A function of the type K <J(9) is a typical example of a corner-singularity 
as arising from the solution of an elliptic boundary value problem. 

Theorem 5.1 (or corollary 5.1) is thus well suited to predict the optimal 
order of convergence (modulo e ) that one can in general expect by dimensional 
reduction of elliptic boundary value problems. 
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A result like this could not have been obtained by using the a-priori 
knowledge of the regularity of solutions to elliptic boundary value problems 
in terms of standard Sobolev spaces. 


i 
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6. A Simple Example of Dimensional Reduction 
Let us consider the boundary value problem 


Au * 0 

in 

]0,l[x3-l,l[ 

u - 0 

for 

x ** 0,1 

3u .. 

an' 8(10 

for 

y - + i 


(n is the outward normal). 

From [133 it follows that the optimal choice of basis functions for dimen- 

o 

sional reduction in this case is the polynomials. denotes the set 

N . Gl 

{ 2, w (x) p. (y) |w.£ H ([0,1])} , where p is a polynomial of degree j . 

j«0 J J J J 


Let u^ denote the projection of u onto V N in the inner product 




1 1 
f 


0-1 


(l± li.aiii) 

\3y 3y 3x 3x / 


dydx 


It is clear that 


inf 

o 


q 


€ V 


N 



2 

H 1 ([0,l3x[-l,l3) 


_< C BCu-Ujj.u-Ujj) £ 


< C 


inf 

o 


q€V. 


N 



2 

H 1 ([ 0, l3x[-l, 1]) 


and hence that the energy error B(u-u^, u-u^) is asymptotically in N equivalent 

1 0 

to the square of the distance (in H ) from u to V N 
If g has the Fourier series 


g(x) - l g, sin kir 
k-1 


x 
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Chen ic immediately follows chat u is given by 


u<x,y) 


00 


l 

k-1 


cosh(kiTy) ^k 
sinh(kir) kir 


sin kirx 


In terms of regularity of u it is not difficult to prove that this formula leads 
to the following three results: 


i) V<*> 0 : l g 2 k 2a < - <b^u€H :3/Z+a ([0,l]x[-l,l]) 

k-1 K 


,3/2+a. 


ii) \/“ > 0: l g l k 2a < - 


k-1 


fue 0 (H 1 ([0,l]x[-l,ll);K R ) fl/ _ 

. RSN ,R>0 6/R> 

I with 0 = a + 1/2 


iii) Va> 0, e > 0: u€ D (^([O.ljxl-l.l]) „ 

R€IN,R>0 0/R » 

with 0 - a + 1/2 + e . 


-A r 2 . 2a < 
^ L 8v k < 


k-1 


We consider two different choices for g 


g(x) - tt/ 4 
g(x) - x(x-l) 


For the first choice of g it follows that 


r 2 20 

Z g v k < » for any 0 < 1/2 , and 
k-1 K 



m oo 






similarly for the second choice 


2 20 

l g. k < « for any 0 < 5/2 , and 


k-1 


2 ,5 


I ^ 


k=l 


Corollary 5.1 together with the regularity results ii) and iii) now 
ensure that 

In the case g(x) « it/ 4 

BCu-Ujj, u-u^) will converge to zero faster than N , y e > 0 , but 

on the other hand slower than N ^ > 0 

In the case g(x) =» x(x-l) 

B(u-u n> u-u n ) will converge to zero faster than N ,ye > 0 , but 

on the other hand slower than N ^ £ e > 0 

Figures 1 and 2 show the actual computed values of B(u-u^, u-u^) as a 
function of N in the two different cases. Note that the asymptotic rate of 
convergence is obtained already for a fairly small number of polynomials. 

For details concerning the computation of the u^'s see [13]. 

Instead of using Corollary 5.1 to obtain information about the rate of 
convergence we could have used the regularity result i) and a 2-dimensional 
version of Theorem 4.1. This way we could at most have predicted convergence 
of the order of N ^ +e and N ^ +€ respectively, i.e. only half the actual 
convergence rate. 

In C133 we considered the same boundary value problem as here, only it was 
on the domain [0,l]x[-h,h] for some h > 0 , and not on [0,l]x[-l,l] . From 


















the computational results there it follows that for a fixed N ^ 2 , 
h h 2 

B(u-u„, u-u ) behaves like h , h -*■ 0 , in the case where g(x) * it/ 4 
" N N 

( u„ is the projection of u onto { £ w . (x) P.(y/h) | w .€ H x ([o,l])} ) . 

N j-0 J J J 

Comparing this to the result obtained here for g(x) = tt/ 4 it is seen that using 

N polynomials, y€C-l,l3, is in some sense equivalent to having a domain of 
2 

thickness 1/N . A similar feature has been noticed by comparison of the 

standard h-version of the F.E.M. with the so-called p-version (cf. [2j). 

In this example we used slight variations of the approximation results 
proved in sections 4 and 5, namely with fixed boundary conditions =0 at 
x * 0,1 . The proofs of these results follow immediately from the proofs of 
the similar results with no boundary conditions. 








In sections 2 and 3, we used some results concerning the eigenvalues and 


eigenfunctions of the boundary value problem 

d d 

- ■%— a -7— u = Xu 
dy dy 

u(-l) = u(l) = 0 

a here is a function in L ([-1,1]) such that 3 a constant a with 

-* o 

0 < a Q <_ a(y) . From the theory of Sturm-Liouville systems it immediately 
follows that the eigenvalues (repeated according to multiplicity) form a sequence: 


0 < X < X, < 
o — 1 — 


< X < X . . 

— m — m+1 


with + « as the only limit point. 

Lemma A.l 

With notation as above 


X t for ® J 4 m' » 

m m 

i.e. the eigenvalues are all simple. 

Proof 

Assume that for some m j* m' X^ " ^ m » • This means that the eigenvalue 
X * X ffl (* X^,) has multiplicity > 1 . Let u and u be two linearly independent 

*Y, 

eigenvectors corresponding to X , and let v ■ cu + du be a nontrivial linear 

dv 

combination with the property that a -j— * 0 for y ■ -1 (such one obviously 




A-2 


exists). The function v is then a solution to the initial-value problem 


d d . 

- 3 — a — v ■ Xv 
dy dy 



0 


in [- 1 , 1 ] 
for y - -1 


and because of the uniqueness of solutions to this problem, it follows that 

v * 0 . Since v is a nontrivial linear combination of u and u this shows 

a. 

that u and u are linearly dependent. We therefore have arrived at a contra¬ 
diction, i.e., X j 4 X , for m j 4 m' 
m m 

□ 

It is well known that there exist constants 0 < and 0 < such that 
C. (nri-1 ) 2 < X < C 0 (m+1 ) 2 . 

By imposing an extra smoothness requirement on a we can obtain a much more 
detailed statement. 


Lemma A.2 

If a £ C 2 ([-l,l]) , then 

X - (ir/i) 2 .(m+1) 2 + 0(1) 

m 1 

f _i /? 

where Z • j (a(y)) dy 
-1 


A proof of Lemma A.2 is found in chapter 4 of [11], and shall not be 
repeated here. 

OO 

Let {u } n denote a sequence of normalized eigenfunctions, u corres- 
m m*»u m 

ponding to X q . Let f^ be given as in section 2, namely 






Lemma A.3 


The function fg has the expansion 


00 



where a ^0 for every m . 
m 


Proof 

That fg has a unique expansion is wellknown. The coefficient 
given by 


1 y 

r r i 

a m " j J I(I) ds U m (y) dy 
-1 -1 


Now assume that for some value of m ■ m a ■ 0 , i.e. 

o m 

o 


i y. 


-i -i 


— 7 —x ds u (y)dy * 0 

a(s) m 0 


Since u (y) 
m 

o 


1 d d 

- — — a 3 — u we get that 

A dy dy m 

m J J 0 

o 


; ; y 

- 1-1 


A) ds [37 * I7 U »J <y) dy ' 0 


Performing an integration by parts this yields 







■ *it, - ' ' ii 




rr™'irJ T» 


A-4 


and the last integral here vanishes due to the fact that u (1) ■ u (-1) » 0 

m m 

o o 

We therefore conclude that 


d 

a — u » u 
dy m m 

o o 


for y * 1 


On the other hand u satisfies the differential equation 

o 


d d 

— a —— u + A u =0 
dy dy m mm 
o o o 


in [-1,1] 


Because of uniqueness of solutions to the initial-value problem, this implies 

that u * 0 . We have thus arrived at a contradiction, meaning that a i 0 
in m 

o 

for every m . 


□ 


Again, by imposing an extra smoothness requirement on a , we obtain a 
very detailed result concerning the decay-properties of the a m ' s 

Lemma A.4 
2 

If a€C ([-1,1]), then ^ constants 0 < and 0 < such that 



for all m 


Proof 

From [ll] p. 176 we get the following asymptotic formula for 


u m (y) 




u m (y) 


- sr «« “• (fcp-«)] 


0((nrt-l)" 2 ) 




where 5(y) ■ j (a(s)) ds and l * £(1) =■ j (a(s)) ds . The function 


-1 


-1 


T is in C and the constants D satisfy 

m 


^ D (independent of m ) such that 


1/D < |D < D for all m 
in — 


Also, O(') here means uniformly in y . Let us now calculate a 


m 


1 y 


m 


-1 


f 0 <y) u m (y)dy - | J ^ds u m (y)dy = 


- 1-1 


.-2s 


\ + x 2 + 0 ( (nrt- 1 ) *) 


I 2 denotes the integral 


m+T °m . | | (a(y ))‘ 1/4 T(()cos (^>1. 5 ) dy . 


By a change of variables from y to 5 and an integration by parts, it 


.-2s 


immediately follows that ^ is 0 ((m+l) ) 


1 ^ is given as 


1 / X 


m 


^ 7 ) ds] (a(y ))‘ 1/4 sin 


-1 \-l 




dy 


By a change of variables from y to £ and an integration by parts we 








IP 


r* 




A-6 


I 


1 



£ (a(l)) 1/4 • 


1 

(m+l)u 


*(-l) m+1 + 0((m+l)' 2 ) . 


This immediately implies the existence of two constants 0 < and 0 < 
such that 




nri-1 


for m sufficiently large. Now combining with Lemma A.3 and possibly changing 
the constants and C2 we get the desired result. 
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